Appearance of Universal Metallic Dispersion in a Doped Mott Insulator 



00 
O 

o 

(N 

<D 

m 



•4— > 

CO 

-4-* 



i 

C 

O 

o 



in 

m 

ON 

O 

00 

o 



S. Sahrakorpi 1 , R.S. Markiewicz 1 , Hsin Lin 1 , M. Lindroos 1 ' 2 , X. J. Zhou 3 - 4 , T. Yoshida 5 , W. 
L. Yang 3 ' 4 , T. Kakeshita 6 , H. Eisaki 3 - 6 , S. Uchida 6 , Seiki Komiya 7 , Yoichi Ando 8 , F. Zhou 9 , 
Z. X. Zhao 9 , T. Sasagawa 3 ' 10 , A. Fujimori 5 , Z. Hussain 4 , Z.-X. Shen 3 - 4 , and A. Bansil 1 
Physics Department, Northeastern University, Boston, Massachusetts 02115 
2 Institute of Physics, Tampere University of Technology, P.O. Box 692, 33101 Tampere, Finland 
3 Dept. of Physics, Applied Physics and Stanford Synchrotron Radiation Laboratory, Stanford University, Stanford, CA 94305 
4 Advanced Light Source, Lawrence Berkeley National Lab, Berkeley, CA 94720 
5 Dept. of Physics, University of Tokyo, Bunkyo—ku, Tokyo 113-0033, Japan 
6 Dept. of Superconductivity, University of Tokyo, Bunkyo—ku, Tokyo 113, Japan 
7 Central Res. Inst, of Electric Power Ind., 2-11-1 Iwato-kita, Komae, Tokyo 201-8511, Japan 
8 Institute of Scientific and Industrial Research, Osaka University, Ibaraki, Osaka 567-0047, Japan 
9 National Lab for Superconductivity, Inst, of Physics, 
Chinese Acad, of Sci., Beijing 100080, China and 
10 MSL, Tokyo Institute of Technology, Kanagawa 226-8503, Japan 
(Dated: September 13, 2008) 

We have investigated the dispersion renormalization Zdisp in I^-xSr^CuCU (LSCO) over the 
wide doping range of x = 0.03 — 0.30, for binding energies extending to several hundred meV's. 
Strong correlation effects conspire in such a way that the system exhibits an LDA-like dispersion 
which essentially 'undresses' (Zdisp — ► 1) as the Mott insulator is approached. Our finding that 
the Mott insulator contains 'nascent' or 'preformed' metallic states with a vanishing spectral weight 
offers a challenge to existing theoretical scenarios for cuprates. 



PACS numbers: 79.60.-i, 71.18.+y, 74.72. Dn 

Under strong electronic correlations the parent com- 
pounds of all cuprate high-temperature superconduc- 
tors assume the so-called Mott-Hubbard insulating state, 
rather than the conventional metallic state. By what 
routes these insulators accomplish the miraculous trans- 
formation into superconductors with the addition of elec- 
trons or holes is a question of intense current inter- 
est, which bears on ongoing debates surrounding the in- 
terplay between electron correlations, magnetism, lat- 
tice effects, and the mechanism of high-temperature 
superconductivity. 1 In this study we consider the clas- 
sic superconductor I^-aSr^CuC^ (LSCO) over the wide 
doping range of x = 0.03 — 0.30, delineating how the 
electronic spectrum evolves with doping for binding en- 
ergies extending to several hundred meV's. Our analy- 
sis indicates that this Mott insulator contains 'nascent' 
or 'preformed' metallic states, which develop finite spec- 
tral weight with doping, but otherwise undergo relatively 
little change in dispersion over a wide doping range. 
Our findings challenge existing theoretical scenarios for 
cuprates. 

We have carried out extensive angle-resolved photoe- 
mission (ARPES) measurements from LSCO single crys- 
tals covering a wide range of dopings, momenta and bind- 
ing energies. Although the incoherent part of the spec- 
trum behaves quite anomalously, we find that many-body 
effects conspire in such a way that insofar as the coherent 
part of the spectrum is concerned, at least its underlying 
dispersion is reasonably described by the conventional 
band-theory picture, significantly broadened lineshapes 
and 'kinks' in the dispersion notwithstanding. Surpris- 
ingly, even with the addition of just a few percent holes in 



the insulator, the full-blown metallic spectrum seemingly 
turns on with little renormalization of the dispersion. In 
particular, the spectrum displays the presence of the tell- 
tale van Hove singularity (VHS) whose location in energy 
and three-dimensionality are in accord with the band the- 
ory predictions. Furthermore, this metallic spectrum is 
'universal' in the sense that it depends weakly on dop- 
ing, in sharp contrast to the common expectation that 
dispersion is renormalized to zero at half-filling. 

The band structure results are based on all electron, 
full-potential computations within local density approx- 
imation (LDA) using the tetragonal lattice structure 2 , 
and include effects of La/Sr substitution within the 
framework of the virtual crystal approximation 3 . The 
ARPES measurements were carried out on Beamline 
10.0.1 at the ALS using Scienta 200, 2002, and R4000 
electron energy analysers for 55 eV light with strong in- 
plane polarization. The energy resolution is 15-20 meV, 
and the angular resolution is 0.3 degrees for the 14 de- 
grees angular mode. All data were taken at T=20K. 

We set the stage for our discussion with the help of 
Fig. 1, which shows typical spectra from LSCO in the 
form of energy distribution curves (EDCs) at two differ- 
ent dopings for a series of momenta. Considering the 
over doped case (upper light green set), we see a coher- 
ent feature dispersing to higher binding energies and be- 
coming broader as one moves away from the Fermi mo- 
mentum hp. This feature sits on top of a substantial 
incoherent background extending to quite high energies 
at all momenta. These basic characteristics are seen to 
persist in the lightly doped sample (lower dark blue set), 
although the greatly reduced spectral weight of the co- 
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FIG. 1: (Color online) Illustrative ARPES spectra as a func- 
tion of binding energy in La2-xSr a: Cu04 for a series of mo- 
menta along the nodal (i.e. V to direction. Results 
from a lightly doped insulating sample (x — 0.03, dark blue 
lines) and an overdoped metallic sample (x = 0.22, light green 
lines) are shown. Coherent spectral peak is seen to disperse 
to higher binding energies as one moves away from the Fermi 
momentum, Icf- 

herent feature in relation to the incoherent part of the 
spectrum is very evident. Our focus in this article is on 
the aforementioned coherent feature in the spectrum of 
LSCO, and especially on delineating the evolution of its 
dispersion with doping. 

The existence of a large, Luttinger-like, metallic 
'nascent' or 'underlying' Fermi surface in LSCO has been 
established in previous studies, culminating in the recent 
systematic analysis of Sahrakorpi et of. 4 and Yoshida et 
al. 5 . In contrast, here we consider spectra over a wide en- 
ergy range of several hundred meV's, show the presence 
of the VHS-a unique feature of the band structure-even 
in the lightly doped insulator, and establish unequivo- 
cally the existence of near-universal metallic dispersion 
in LSCO. These 'nascent' Fermi surfaces and dispersions 
are well defined despite the difficulties of identifying fea- 
tures in the face of loss of spectral weight as the pseudo- 
gap develops with underdoping. We emphasize that our 
focus is on what we may call the 'gross' spectrum. In 
other words, we are not concerned with the fine struc- 
ture in the electronic spectrum associated with the well- 
known low energy kinks 6 , the recently discovered fea- 
tures at higher energy scales, 7-11 or superconducting 12 , 
or other 13 leading-edge gaps. 

The top row of Fig. 2 shows ARPES intensity maps 
from a lightly doped sample of LSCO (x = 0.03) for a se- 
ries of binding energies. Cross-sections of the correspond- 
ing constant energy (CE) surfaces in the {k x ,k y ) plane 
computed from the band structure of LSCO, superposed 
at k z = (magenta lines) and k z — 2n/c (black lines) 
indicate the expected broadening of the ARPES spectra 
associated with interlayer coupling. At zero binding en- 
ergy, i.e. the Fermi energy Ep, such CE contours give 



the projection of the 3D FS of LSCO on to the (k x , k y ) 
plane. Notably, the momentum region enclosed by these 
CE contours defines the region of allowed ARPES tran- 
sitions, modulated by the effect of the ARPES matrix 
element. 4 ' 14,15 At low binding energies, the CE surface is 
seen from panels (a-c) to be hole-like around the X(ir, it) 
point for all k z values. In contrast, at high energy in 
panel (e), after the VHS has been crossed, the CE sur- 
face becomes completely electron- like centered at T. The 
transition from a hole- to electron-like CE surface does 
not take place abruptly because the VHS possess a signifi- 
cant 3D character, extending from 85-140 meV in binding 
energy. 

The evolution of the experimental ARPES intensity 
pattern with binding energy in the top row of Fig. 2 
clearly follows that of the projected CE surfaces. In par- 
ticular, the spectral intensity remains confined mainly 
within the boundaries of these projections as expected, 
and with increasing binding energy, the intensity first 
spreads towards the M(n, 0)-points and then moves away 
from the M-points along a perpendicular direction, very 
much the way the CE surfaces transition from being hole- 
to electron-like. Moreover, first principles ARPES com- 
putations show that under the combined effects of the 
matrix element and k z dispersion, the spectral inten- 
sity develops the characteristic 'wing-like' shape seen in 
Fig. 2, and that the spectral weight grows rapidly in the 
antinodal region as the VHS is approached. 4 These re- 
sults leave no doubt that metallic states, including the 
presence of the 3D VHS, appear in the spectrum of the 
insulator with the addition of only a few percent holes, 
even though spectral broadening and incoherent back- 
ground make it hard to see this directly in the individ- 
ual EDCs. The observed location in energy and three- 
dimensionality of the VHS is well-described by the con- 
ventional band theory picture, indicating that the ener- 
gies of these metallic states undergo little renormalization 
in the lightly doped insulator. The evolution of the VHS 
over a denser binding energy mesh can be seen further in 
the three movies for x = 0.03, 0.07, and 0.12 provided as 
supplementary material. 16 The characteristic dispersion 
of the VHS is quite recognizable in the emission maps of 
the top row of Fig. 2. and in the supplementary movies. 

We demonstrate next that the aforementioned metallic 
dispersion is only weakly dependent on doping. Note that 
if this is true then the main difference in going from one 
doping to another would be a shift in the Fermi energy 
needed to accommodate the right number of holes in the 
filled portion of the band structure. That is, topologies 
of the CE surfaces and the associated emission spectra 
would be comparable for various dopings except for a 
rigid shift of the energy scales. That this is indeed the 
case is shown by the results for emission from the Ep for 
x = 0.05 — 0.30, presented in the bottom row of Fig. 2. 
For example, based on parameter free LDA computations 
for x — 0.03 and x = 0.30, the Ep for x — 0.30 is lower by 
190 meV than for x = 0.03, so that the ARPES map for 
emission at a binding energy of 190 meV from x = 0.03 in 
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FIG. 2: (Color) Experimental ARPES intensity maps in LSCO are compared with the corresponding cuts in the (k x , k y ) plane 
through the theoretical constant energy (CE) surfaces for k z = (magenta lines), and k z = 2n/c (black lines). Top Row: 
x = 0.03 with binding energy varying from 22 meV to 190 meV. Bottom Row: ARPES maps for emission from the Fermi 
energy for dopings varying from x = 0.05 to x — 0.30. 



panel (e) can be compared with that for emission from 
the Ep for x = 0.30 in (j). In this vein, the binding 
energies in the various panels of the top row of Fig. 2 
for x — 0.03 have been chosen to match the Ep shifts in- 
volved at the doping levels considered in the panels of the 
bottom row. Good accord is seen in all cases. We empha- 
size that the energy shifts corresponding to the respective 
panels of the top and bottom rows in Fig. 2 are based on 
parameter free first principles computations. These and 
other similar comparisons among spectra taken at differ- 
ent binding energies and doping levels show clearly that 
LSCO is characterized by a near-universal metallic dis- 
persion despite dramatic changes in the lincshape due to 
interactions. This is also true for the theoretical disper- 
sions, although slight doping dependencies can be seen 
for example by comparing the CE surfaces in the top 
and bottom frames (c) and (h) in Fig. 2. 

Further insight is provided by Fig. 3, which shows plots 
of spectral intensity as a function of binding energy along 
the antinodal line (top row) and the nodal line (bottom 



row) for six different dopings. The corresponding en- 
ergy bands at k z = 0, tt/c and 2ir/c are overlaid in each 
panel. These bands are seen to differ substantially for 
different k z values along the antinodal line in the upper 
panels but are virtually indistinguishable in the lower 
panels along the nodal line. The VHS, which is marked 
by the extremum of the band along the antinodal line, 
is spread over 40 meV at x = 0.03 doping due to the 
effect of k z dispersion, and its center lies at 110 meV be- 
low Ep for x — 0.03, but moves to 110 meV above the 
Ep for x = 0.30 with a spread of 70 meV. Along the 
nodal direction, the theoretical band follows the exper- 
imental peak positions given by the black dots reason- 
ably well, although the experimental points are shifted 
to the right compared to theory in most cases, indicat- 
ing a slight deviation of FS shape from LDA. Along the 
antinodal direction, however, the spectral peaks are sub- 
stantially broader due to k z dispersion and also possible 
many body interactions, although most of the spectral 
weight lies within the k 2 dispersed bands as expected in 
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FIG. 3: (Color) ARPES intensity maps along the antinodal (top row) and nodal (bottom row) directions over the doping range 
x — 0.30 — 0.03. The corresponding computed band structures are also plotted for three different values of k z : k z = (magenta 
solid line), k z — ir/c (magenta dashed line), and k z — 2n/c (black solid line). Black dots mark positions of the peaks in the 
experimental spectra. Note that nodal data in the lower row are plotted on an expanded horizontal scale in order to highlight 
relatively small differences between gross theoretical and experimental dispersions. For each doping, the nodal and antinodal 
intensities are normalized to a common maximum. 



Doping x = 0.03 0.07 0.12 0.15 0.22 0.30 

Z dlsp (nodal) 1.0 0.8 0.6 0.6 0.6 0.6 

(1.2) (1.2) (0.8) (0.8) (0.7) (0.6) 

Z dlap (antinodal) 1.1±0.2 1.3+0.4 >0.6* >0.6* >0.7* >0.7* 



TABLE I: Estimates of dispersion renormalization factors 
Zdisp in LSCO for different dopings x in relation to the LDA 
values. Nodal Zdisp values in parenthesis are an estimate of 
the upper limit. Stars in the second row denote that for these 
dopings the values refer to the region in the vicinity of the 
antinodal point as discussed in the text. 



a quasi-2D system. 4 ' 14,15 

We have examined the renormalization of dispersion in 
relation to the LDA values along the nodal and antinodal 
directions as a function of doping, and thus obtained the 



associated renormalization factors Zdisp- These results 
are summarized in Table I. Our study provides new in- 
sight into the value of Zdisp along the antinodal direction 
because the VHS is a very robust feature of the LDA 
band structure. By determining the position of the VHS 
in the experimental spectrum, and comparing this posi- 
tion with that expected from the LDA, we can uniquely 
determine the overall renormalization of the spectrum 
along the antinodal direction with respect to the LDA. 
At x = 0.03 and x = 0.07, the VHS lies below the Fermi 
energy for all k z values, and in these two cases, by analyz- 
ing the supplementary movies 16 using different renormal- 
ization factors to scale the LDA bands, we have obtained 
Zdisp (antinodal) values of ~ 1 within the uncertainty 
shown in Table I, so that the LDA bands are essentially 
unrenormalized. For x = 0.12 and higher dopings, part 
or all of the VHS lies above the Fermi energy, so that we 
were only able to estimate a lower limit for Zdisp in the 
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FIG. 4: (Color online) Reproduction of Fig. 1 in Zhou et al. 17 
Dark blue line drawn on the x = 0.03 dataset illustrates how 
the nodal Zdisp value is defined here, while the light green line 
gives the high energy slope used to estimate the corresponding 
upper limit on Zdisp (nodal) as discussed in the text. 



vicinity of the antinodal point. There are of course no 
filled states at the antinodal point once the VHS moves 
above the Fermi energy. 

Turning to the nodal direction, we note that here the 
effects of k z dispersion are small, and our analysis is con- 
sistent with results available in the literature. For com- 
pleteness, however, we have estimated the gross Zdisp 
(nodal) values from the nodal spectra given in Fig. 1 of 
Zhou et al. 17 , reproduced here as Fig. 4. In order to 
gain a handle on the 'gross' or underlying dispersion ex- 
clusive of the low energy kink (see arrow at 70 meV in 
Fig. 4), values of Zdisp (nodal) given in Table I are ob- 
tained from the slopes of the straight lines joining the 
point at the Fermi energy with that at 200 meV binding 
energy, i.e. by lines such as the dark blue line drawn for 
the x = 0.03 case in Fig. 4, and comparing this slope to 
the corresponding LDA values. Since the LDA values of 
the nodal radii do not exactly match the measured val- 
ues, we have estimated Zdisp (nodal) by comparing LDA 
and experimental slopes. We also show in parentheses 
the values of Zdisp (nodal) obtained from the slopes of 
the high energy part of the spectrum, as given for exam- 
ple by the light green line for x = 0.03 in Fig. 4, which 
provides an estimate for the upper limit of Zdisp (nodal) . 

Interestingly, at the lowest doping of x = 0.03, the 
nodal as well as the antinodal renormalization factor is 
seen to be ~1. In the optimally and overdoped regimes, 
the nodal renormalization factor is ^0.6, while the value 
of the antinodal factor is estimated to be greater than 0.6. 
In the underdoped x — 0.03 case, the renormalization of 
states in the antinodal and nodal directions is roughly 
similar, but this is less clear at higher dopings. These 



results are surprising since we might have expected the 
LDA to provide a reasonable description in the overdoped 
regime, and to be renormalized greatly in the underdoped 
case. 

Even though we have shown that the gross dispersion 
up to several hundred meV's is described quite well by 
the band theory picture, we emphasize that this does 
not mean that the spectrum of LSCO is conventional 
in nature. As shown in Fig. 1, the coherent spectral 
weight of these dispersive features fades away for under- 
doping where the pseudogap and polaronic effects kick 
in. The spectral weight is transferred to an incoher- 
ent feature 18-20 , reflecting the importance of many-body 
physics. Our analysis thus suggests that the spectrum 
of the insulator already contains 'preformed' or 'nascent' 
metallic states, which possess zero spectral weight in the 
half-filled case. With doping, these states develop finite 
spectral weight, but otherwise undergo relatively little 
change in their basic character. 

It is interesting to study the relationship between the 
renormalization of the dispersion and spectral weight 
within the simple Green function formalism. The Green 
function can be approximately written as G = Z u /[uj — 
Z u e k /Z k ], with Z w = 1/(1 - dE/du), and Z k = 1/(1 + 
dYj/dek) 1 where S is the self-energy. While the spectral 
weight renormalization is associated with Z u , the band 
dispersion renormalization is given by Zdisp = Z^/ Z k . It 
is striking that the spectral weight in LSCO is renormal- 
ized greatly but the dispersion is not. This is not consis- 
tent with a momentum independent self-energy, since the 
renormalization of spectral weight would then be accom- 
panied by a similar renormalization of dispersion, sug- 
gesting that the self-energy possesses momentum as well 
as energy dependencies such that the effects of the asso- 
ciated renormalization factors, Z^ and Zk, 21 conspire to 
approximately cancel each other in renormalizing disper- 
sion. This picture however implies that the heat capacity 
would be renormalized only weakly with doping, which 
is not the case experimentally. 22 

Insofar as theoretical models of strongly correlated sys- 
tems are concerned, a focus of attention in the literature 
has been the t — J model, which has been used to de- 
scribe the physics of the Hubbard Hamiltonian in the 
large U regime. One class of t — J models finds gapped 
insulating solutions at half-filling, which evolve with dop- 
ing to yield small Fermi surface pockets 23 . This route 
appears to be followed in Nd2- x Ce x CuOi for electron 
doping 26 . But the present results do not show clear 
evidence for small pockets or other strong violations of 
Luttinger's theorem 27 . On the other hand, metallic solu- 
tions with a large Luttinger-like Fermi surface have also 
been reported 12,28 , where the spectral weight Z^ is renor- 
malized to zero at half-filling. These solutions are closer 
to our observations, even though the band width is found 
to be strongly renormalized near half-filling, Zdisp ~ J/t- 
A recent variant of the t — J model succeeds in captur- 
ing both the metallic and insulating aspects of the spec- 
trum simultaneously 29,30 . However, to our knowledge, 
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our finding of a large Fermi surface and a large disper- 
sion, which is weakly dependent on doping, does not fit 
well within the framework of the currently available sce- 
narios based on simple Hubbard or t — J models or naive 
application of the classical polaron theory. 

In conclusion, we have shown that metallic states ap- 
pear in the electronic spectrum of high temperature su- 
perconductor L^-zSr^CuCU at all doping levels, rang- 
ing from the lightly doped insulating to the overdoped 
metallic regime. The gross dispersion of these states is 
not only near-universal in that it is weakly dependent 
on doping, but it is also conventional in that it is rea- 
sonably described by the standard band theory picture 
up to binding energies of several hundred meV's. How- 
ever, the metallic states suffer substantial loss of spec- 
tral weight with underdoping and in this respect behave 
quite unconventionally. Our findings challenge existing 



theoretical models of the cuprates and indicate the com- 
plexity with which many-body physics plays out in this 
fascinating material. 
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